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Summary. — In the 90s of the last century high statistics experiments with fully 
equipped 4n detectors have lead to a better insight in the spectrum of hadrons. 
In particular the finding of crypto-exotic and exotic states tremendously im- 
proved the experimental situation in meson spectroscopy. All this was possible only 
with sophisticated analysis methods like the decomposition of measured phase-space 
distribution into partial waves and to express the partial waves in terms of compli- 
cated dynamical functions. This paper gives an introduction about the concepts 
and formalisms involved. 

PACS 11.80 - Partial wave analysis. 



1. Introduction 

l'l. Goals. - The spectrum of hadrons, mesons and baryons is the result of bound 
states of quarks like qq and qqq respectively. Although the quark-model was moti- 
vated partly by the structure of the spectra and the symmetry and decay properties 
of the hadrons, the gauge field theory of strong interactions, Quantum Chromodynamics 
(QCD), has a very strong coupling constant at very low momentum transfer (e.g. low 
masses) so that there are no firm predictions in the classical field of light quark spec- 
troscopy and medium energy physics. Therefore there is no possibility to map out the 
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states from first principles (if we leave lattice gauge theory out for the moment). To 
understand the spectrum of light hadrons it is important to collect the properties to de- 
rive effective models to be used in that field, which may tie up to the very high energies 
where perturbative QCD works well. In order to uncover the spectrum it is necessary 
to investigate the scattering and the decay of hadrons and to identify the intermediate 
states throughout the whole reaction. 

Lattice gauge theory could be a way out of the problem, but since the results reflect 
a measurement on the lattice (with very limited precision and a lot of assumptions so 
far) it is not suited to actually understand the principles and the physics behind the 
measurements, e.g. an actual state and/or it's properties. 

The basic task is to find all resonances, with their static properties like mass, width, 
spin and parities. Since SUf(3) [1^ provides a lot of relations among the decays of pure 
states (without Fock states) without complicated final state effects it is important to 
measure the partial decays widths as well. This is a very demanding task, since a lot 
of resonances overlap. In addition complicated production processes or scattering with 
many waves in the intermediate state complicate the situation. To disentangle the waves 
and to identify resonances and their actual yield. 

1'2. Technique. - The techniques are manifold and we should to concentrate only on 
two basic approaches. The first one is scattering of hadrons, usually on a proton or a 
deuterium target, depending on the desired process. Typical examples are 

• 7rN scattering with or without charge exchange (GAMS at CERN, E852 at AGS), 

• 7 N scattering (CEBAF, MAMI, Elsa, Graal), 

• 7rp or pp in the central region (WA76, WA92, WA102 at CERN, E690 at FNAL), 

• pp near meson production thresholds (WASA at Celsius, Anke and TOF at COSY), 

• pp in flight (Crystal barrel, Jetset and PS185 at CERN, Panda at GSI). 

These reactions involve usually a lot of partial waves. Single or double polarization ex- 
periments, low excitation energies and/or a selection of specific exclusive channels which 
impose reasonable constraints on the reaction are required to enable an unambiguous 
decomposition of the system. 

The second one is where the initial system is at rest and/or the remaining part of the 
reaction is not of interest like 

• pN reactions at rest into many body final states (Asterix, Crystal barrel and Obelix 
at LEAR), 

• K° and Jf± decays (NA48 at CERN, Kloe at Da^ne, kTev at FNAL), 

• 0(1020) decays (Kloe at Da^ne, VEPP at Novosibirsk), 

• D and D s decays in high energy reactions (photo-production at FNAL, Babar/Belle 
at PEP-2/KEK-B, CLEO-c at CESR), 
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• J/V> decays (Marklll at SLAC, DM2, CLEO-c at CESR, BES at BEPC). 

For all experiments of this type it is important that spin densities are well known. Oth- 
erwise a full decomposition is not possible. In that case additional assumptions have to 
made to relate the different parameters. 

13. Methods. - 

1'3.1. The Partial Wave Approach. 



Motivation. Partial waves are easily introduced in a scattering process. To get a first 
impression we start with Schrodingcr's equation 



(1) 
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with k = — = u— and the reduced mass a = . The incident wave can be 

n n mi + mi 

expressed like ^>i(r, 1?, ip) = e tkz and we assume a vanishing potential V(r) = 0. Then 

we can expand the initial state \i) in terms of Legendre polynomials Pi, thus separating 

angular and radial wave function 



(2) 



= =J2Ui(r)Pi(cos#). 



The scattering wave function ^5 is the difference between incoming and outgoing wave. 
We parametrize Ui in terms of a phase Si and an inelasticity r\i, which are motivated 
from our knowledge of resonance curves, where the phase moves from — -| to ^ and the 
inelasticity carries all dissipated probability (e.g. into other channels): 
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The total cross section can then be written 
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the so called T-matrix. We will discuss it's properties later in more detail (see sec. 13'2.4^1 . 
To visualize the properties of T Argand plots are quite useful, since they provide a check 
on the unitarity of T (which ensures that the probability of a reaction does not exceed 
unity). The Argand plot is a graph of the imaginary part of T versus the real part of 
T (see fig.^l. It is seen easily from this plot that the phase Si is not the same as the 
angle <p which apears in the polar representation of complex numbers. As an example, 



flm(T) 




I = \|A|/=I 



Re(T) 



Fig. 1. - Argand plot definitions. 

hg-El shows the behavior of T for a relativistic Breit-Wigner function. If the inelasticity 
is zero, e.g. there is only one elastic channel, the Argand circle has a constant radius of 
2 and crosses the imaginary axis at one. From eq.0]we see that the angular amplitudes 
(Pi) and the dynamic amplitude (Ti) factorize. 

1*3.2. The Isobar Model. The last preparative step to write down an amplitude is 
to make assumptions how particles are grouped to construct a decay /reaction chain. An 
empirical approach is the isobar model. It assumes that all subsequent decays appear 
to be two-body reactions. This model seems to work extremely well in very different 
environments and for most hadrons (exceptions may be lo — > ir + ir~ir° and r/ — > 37r). The 
type of valid reactions is sketched in fig.[2|a, while fig.|Hb shows a more complicated reaction 
involving rescattering. The isobar approach will not work easily in those circumstances. 
In some cases it can be retained if the rescattering process is refactorizable in terms of 
many isobar reactions. This involves usually very many parameters and may not lead to 
a sufficient description. In those cases model assumptions have to be made. In the isobar 
model, the two-body decay into particle 1 and 2 factorizes (fig. El completely from the 
recoil system 3, which might decay as well. Any node in the decay tree is represented by 
the same isobar-like amplitude. The recoil systems are only involved in maintaining the 
conservation of angular momentum and spin-projections, which imply summation over 
unobservable quantities. 

1*3.3. Construction of the Amplitude. The full amplitude for each node of the decay 
tree consists of a dynamical part XJ and an angular part which we will call Ri. If we 
deal with strong or electromagnetic interactions also isospin is conserved, which means 
that we have to include an isospin related part so that decays with different isospins 
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Fig. 2. - Simple relativistic Breit-Wigner. (a) Argand, (b) phase 8 and (c) intensity plot. 



I-/T3) can be related to each other. This is important if the same charged particle occur 
with different charge in various parts of the decay tree. With these ingredients, the full 
amplitude is defined via 

(6) /(j,/ 3j s > n) = /j(/ > /3)T,(«)i2,(n) 

One should note, that after combining the amplitude for all nodes in the decay tree, 
all conservation laws (like J, /, and their respective projections) have to be taken into 
account, which results in a sum over all unobservable values. 




Fig. 3. - Production process with propagation (a) and rescattering (b). (a) can be treated easily 
with an isobar method, while (b) can not. 
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Fig. 4. - Isobar definition in a chain of subsequent decays. 9 defines the angle for angular decay 
distributions. 



Example: Isospin Relations in pp(0 + ,1 ) — * pn. The initial pp state has either 
/ G (J FC ') = l-(0-+) (called ^o) or I ^) = 0~(1— ) (called 3 Si). The two final 
state gammas have I G (J PC ) = + (l ). Both initial states can decay to pw but this 
case illustrates how important isospin Clebsch-Gordan coefficients are, since 

p°tt° -> (ioioioo) = ~y| 

P°tt q -> (1010|10) = 

(7) pVF (1(±1) 1( T 1)|Q0) = 

p ± ttT _> (i(±l) 1( T 1)|10) = ± 



3 

2 



it is evident that 1 Sq has destructive p ± /p°7r° interference, while for 3 Si only 
contributes and p°7r° does not exist. 



2. Spin formalisms 

2'1. Preface. - There are various spin formalisms available. In principle there are 
three basic types: 

• tensor formalisms, in non-relativistic (Zemach) or covariant form, 

• spin-projection formalisms, where a quantization axis is chosen and proper rotations 
are used to define a two-body decay 

• formalisms based on Lorentz invariants (Rarita-Schwinger) where each operator is 
constructed from Mandelstam variables only. 

The tensor formalisms are very fast algorithms if waves with small angular momentum 
are involved but is getting very complicated if higher waves are present and/or a lot of 
subsequent decays occur. Using the Lorentz invariants is not really a formalism, since one 
has to construct an amplitude with the proper transformation and symmetry properties 
as the actual particle and wave to be represented. As for the tensor formalism this might 
be simple and extremely elegant for waves with low angular momentum but is virtually 
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impossible for a complicated decay cascade. The Zemach formalism as a representative 
of tensor formalisms is discussed in sec. 12'21 

If we restrict ourselves now to the other formalisms the key steps in the specification 
of a scattering formalism are as follows: 

• the definition of single particle states of given momentum and spin component 
(p-states) , 

• the definition of two-particle ^states in the s-channel center-of-mass system and 
of amplitudes between them, 

• transformation to states and amplitudes of given total angular momentum In- 
states) , 

• symmetry restrictions on the amplitudes, 

• formulae for observable quantities, 

• specification of kinematic constraints. 

The basic three formalisms are known as helicity, transversity and canonical (historically 
known as orbital) formalisms. Their basic properties are summarized in tbl.[I] Since they 
have different symmetry relations they can be used in different fields. While the helicity 
formalism is applicable to most spectroscopy experiment, it is the transversity formalism 
which conserves parity and can therefore be applied for partial wave decomposition in CP 
measurements. The three formalisms differ mainly in the choice of the spin quantization 





Helicity 


Transversity Canonical (Orbital) 


property 




possible / simplicity 


partial wave expansion 


simple 


complicated complicated 


parity conservation 


no 


yes yes 


crossing relation 


no 


good bad 


specification of kinematic constraints 


no 


yes yes 



Table I. - Comparison of various spin formalisms. Depending on the choice of the quantization 
axis, the amplitudes have different symmetry properties. 



direction for each particle. 

In the helicity formalism each particle spin is quantized parallel to its own direction 
of motion so that its projection, the helicity A is diagonal. 

In the transversity approach, the component r normal to the scattering plane is used. 

In the orbital (or canonical approach) the component m in the incident z-direction is 
diagonal. 

The merits of each approach arise largely from the invariance of the projections. These 
properties are made precise by the definition of the particle states. It is convenient to 
define them by explicit Lorentz transformation from the spin states m of a particle at rest. 
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Detailed discussions are found in the original papers O El ID • The essential principle in 
all cases is, that the spin projection of a particle is defined in its rest frame. The general 
single particle state is then defined by applying a definite sequence of boosts L(jp x ,p y ,p z ) 
and rotation R(a, (3, 7) to the system. We define 



(8) 


*A = 


|P,A> 


= %, 0,-^)5(0,0, p)\m) 


= H(p)\X) 


(9) 


*r = 


\P, T ) 


= IpA)A a ' t = AH(p)A- 

A 


- 1 \r)=f\r) 


(10) 


*m = 


\pm) 




= R- l (<j>,e,-4>)H{$)\m) = 6\m 



X 



It is clear that it is necessary to calculate the effect of a general Lorentz transformation 
on the states \p£). While hp = p' the corresponding state is not since the states 

are being always obtained from the rest state by a definite sequence X(p) and these do 
not in general commute. However, since both have the same value for the momentum 
in a particular system they can both be transformed back to the rest frame by X~ 1 can 
there they can differ at most by a rotation 

(11) L\pO = X(Lp)R(L,p)\00 = Y, \Lp,e)D J ei (r) 

with the Wigner rotation given by 

(12) R = X- 1 (Lp)LX{p). 

The two-particle states are then defined essentially as direct products of single particle 
states, but phase space factors may be included to make the form in the center-of-mass 
system convenient for various purposes. More details are discussed in the particular 
sections about the different formalisms (see sec. l2'3.l|) . 

2'2. Zemach Formalism. - 

2'2.1. Formalism. The Zemach formalism was originally developed for the investiga- 
tion of the K° decays into 37r [5]. The basic concept is that every angular momentum 
involved in the reaction is represented by a symmetric and traceless tensor of rank ^in 
3-dim. phase-space. The tensors A for spins up to two are 

1 = A° = l 

(13) 1 = 1 A 1 (q)=q 



1 = 2 A*{q) = \ 
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with 



/ <?i \ / qiPi qiP2 qiP3 

(14) 9*-^= <?2 ( Pi P2 P3 ) = Q2P1 (12.P2 (l2PZ 

\ Q3 J \ <?3Pl 93P2 <?3P3 

or with all indices 





I = 


A° = 


= i 


(15) 


I = 1 


A\ = 


= q t 




I = 2 


4 : 


3 

~ 2 



The coupling of spins and/or angular momenta, like spin s of a particle and angular 
momentum I relative to another spinless particle is then done by multiplying the tensors 
and contraction of the resulting tensor. 

2'2.2. Example: pp(0 h ) — * /2vr° — > tt + tt^tt°. The two step process pp — > /27r° 
and /2 — ► 7r + 7r~ may serve as an example to illustrate the method. For p being the 
momentum between 7r° and and q the momentum in the subsequent fi — > 7r + 7r~ 
two-body decay the construction of the amplitude leads to 



(16) 



^4° _ Af^o^A^+^-^t ^ikSji^ 

unpolarizcd 

= X] A 2 f 2 na 1 ijA^.+ n - jkl 5ik5ji 

i,3 



with 



(17) 
(18) 



A 2 



\ViPi ~ \\Pi?$v 

3 1, ,2; 



Combining eg. 1181 and eq.El w ith eg .1161 gives 













2^1% J 







= ~(<f • p) 2 - \tf - \tf + 3jkfbT 

(19) =\{q-p) 2 -\ff 
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and finally the angular distribution is calculated by the squared amplitude 



(gpcostf) 2 - -q 2 p 2 



(20) 



J(cos^-i) 2 = P 2 °W 2 . 



2'3. Canonical and Helicity Formalism. 
2*3.1. Few Particle States. 



Canonical (Orbital) Formalism. 



Single Particle States May \jm) be an at-rest state with spin j and spin projection 
m on to the z-axis in an euclidean system with coordinates (x, y, z). The single particle 
state with a momentum p is constructed through a pure Lorentz transformation Lp 



(21) 



\P, jm) d = Lp\jm) 



R L z pR 1 \jm) 



where Rq = Ro{tp, 0) rotates the z-axis (like in fig-EJi) in the direction of the momentum 
p (ep — Ro(f, tht 7 0)e z . In the first step the momentum vector is rotated via R^ 1 in the 



(a) 




(b) 




Fig. 5. - Definition of angles in canonical (a) and helicity formalism (b). 



z-direction. Secondly the absolute value of the momentum is Lorentz transformed along 
z and finally the z-axis is rotated to the p direction via Rq. 

If rotation of the single particle state is derived from the properties of the rotation 
group. One obtains 

(22) %m)=^I)4,j£p;m) 

ra' 

where D J mlrn (R) are the Wigner D-functions for the rotation R. From eq.|23we see that 
the canonical states transform like at-rest states \jm). 
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Two- Particle States Two-particle systems 1+2 in a rest-system J with respective 
spins can be constructed directly with the help of single particle states. One particle s 
has the momentum p s and the other partilce t has the opposite momentum —p s . Apart 
from that the at-rest system of the two particles depends only on the stereo angle Q s 
and the direction of p s . Using cci.l21lone gets 



(23) 



srrietmt) = k 



L J)s^ \sm s )L Jh^ \tm t ) 
(E s ,p s ) (E t ,-p s ) 

The normalization k can be derived from the single particle states: 

i fpT i 



(24) 



47T V tnj 47T 



where mj is the invariant mass of the state J and p s the invariant phase-space factor. 

To get a practical recipe it is necessary to couple the angular momentum L and the 
total spin S = s + t. Thus the first step is the coupling of the single particle spins to £ 



(25) 



Sms) = /J (sm s tm t \Sms)\fl, sm s tmt) 



where (sm s tm t \Sms) is a Clebsch-Gordan coefficient. A state with a particular total 
angular momentum L is then 

(26) \Lm L Sm s ) = J dnY„ L lL (Q)\n, Sm s ) 

Finally we couple L and S to J: 

\JMLS) = ( Lm L Sm s \JM)\Lm L Sm s ) 

m L ,ms 

(27) = (Lm L Sm s \JM)( S m s tm t \Sm s ) J nYi L {n)\n s , S m s tm t ) 

ttil ,ms ,m s .mt 

with Y^ L (fl) being the spherical harmonics. With these definitions we get the following 
completeness relation 



(28) 



1= Y \JMLS)(JMLS\. 

J,M,L,S 



The normalization of the canonical states is then given by 

(29) (n s > , s'm s 't'mt'\n s , sm s tm t ) = 5(fl s > - fl s )o'ss>5tt>o~m s m s ,S mtmt , 

(30) (J'M'L'S'IJMLS) = 5jj,8 MM ,5 w 8 ss ,. 
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Helicity Formalism. 

Single Particle States The construction of states is similar to the procedure for 
canonical states. May \jX) be an at-rest state with spin j and spin projection A on to the 
z-axis in an euclidean system with coordinates (x,y,z). The single particle state with a 
momentum p is constructed through a pure Lorentz transformation Lp plus a rotation. 
This is due to the fact, that the z-axis is firstly rotated to the direction of p. Secondly the 
result is Lorentz transformed. Therefore the new z-axis, z', is parallel to p (see fig -El) ) • 

(31) \p,j\) d H Lp\j\) 

(32) = R L zP R^\jX). 

The second line of ea. 13 II follows from ea .12 II and the unitarity of rotations. If the helicity 
state is rotated, then the momentum p is rotated, but the helicity is invariant, since the 
quantization axis is rotated as well. Therefore 

(33) R\p,X) = \Rp,X). 

Since eq.[23 and 123 are complete, it is possible to represent the helicity state in the 

r*Fmrmir*al ViasiQ pmH virp vprsa- 



canonical basis and vice versa 

(34) \p,X)=J2 D n^ R 0M 



m) 



Two-Particle States Analogue to the construction of the canonical two-particle state 
ea.l?Tl the helicity equivalent is constructed using the states ea.l3*Tl 

def ^ 

\tt s ,sX s tX t ) = KR [L z p s \sX s )L z p t \tX t )] 

(35) = R (Q s )\n= (0,0), sX s tX t ) 

Using a similar procedure to couple all spins one obtains 

(36) \JMX s X t )=Nj J dn Djl Xg _ Xt \n,sX s tX t ) 
with the normalization factor 



The choice of Nj was made, so that we get a simple completeness relation 



(38) 



1= \JMX s X t )(JMX s X t \ 

J,M,\ s ,\ t 
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The normalization of the helicity states is then given by 

(39) (n s ,,s'\ s ,t'\ t ,\ci s , s x s tx t ) = s(n s , - n s )d ss ,s tt ,s XsK ,d XtXt , 

(40) (J'M'\ s ,\ t ,\JM\ s \ t ) = 6jj,6 M M>8x 3 x s ,Sx t x t , 

2'3.2. Decay Amplitudes. The two-particle states can now be used to derive two-body 
decay formulae in the two formalisms. 

Canonical (Orbital) Formalism. May J be the at-rest system of the decaying state | JM) 
with arbitrary, but well defined spin J and projection M . The decay amplitudes is derived 
by using the two-particle states and summing over all unobservable spin-projections 

(41) A= (p s ,sm s \(-p s ,tm t \M\JM) 

with M. being the unknown decay operator. To get to the final formula we need the 
following relation 

(tt s ,sm s tm t \JMLS) = ^ (Lm L Sm s \JM)(s'm s > t'm t >\Sm s ) 

niL ,ms ,m s * ,m t * 



(42) J dn s ,Y^ L {n s ,)(n siS m s tm t \n s ,,s'm s ,t'm t ,} 

= Y {Lm L Sm s \JM)(s'm s < t'm t >\Sms) 

niL ,ms ,m s i ,m t i 

(43) = (Lm L Sm s \JM){sm s tm t \Sm s )Y^ L {Q s ) 

m L ,m L 

Now we include the two-particle states and obtain 

47T 



(44) A J m ™ mt = —(n s ,sm s tm t \M\JM) 



= y^(n s ,sm s tm t \JMLS)^={JMLS\M\JM) 



W V" y/Iira'l s (n s7 sm s tmt\JMLS) 



def 

L,S 

def 



(45) =' J2 ^a J LS (Lm L Sms\JM)(sm s tm t \Sm s )Y^ L (n s ) 

L.S,niL,ins 

with the canonical partial decay amplitudes a J LS 

(46) a J LS d = f X [^-UMLS\M\JM) 

V Ps 



14 



K. Peters 



which contain the physical matrix element of operator M. 

Helicity Formalism. The decay helicity amplitude is derived in the same way as eg. 1411 

(47) A = (p s ,sK\(-Ps,t\ t \M\JM). 

Using 

(n s ,s\ s t\t\JM\ s ,\ t ,) = Nj J dn' D J A i Xsi _ Ki (n')(n s , sx s t\ t \fi s , , s ' \ S 't' \ t ,) 
= Njj dn s >DfZ !Xsr _ K ,(n S ')5(Q, s > -n s )6 ss ,6 w 5 Xll \ a , 

(48) = NjD{l K ,_ Xt ,(n s ) 
and inserting the two-particle states we get 

(49) A{\ = — {Q s ,sX s tX t \M\JM) 

Ps 

4:7r 

= 2J (Cl a) sX s tXt\JMX ar mt>) — (JMX s >Tnt>\M\JM) 
= J—(2J + \){JMX s m t \M\JM)Di; x ,_ x (n.) 

V PS .st 

(50) =A r j/A s A t ^;, Aa ,-A t ,(^) 
with the helicity amplitudes 



(51) Njf XsXt = Jy(2J +l)(JMX s m t \M\JM) 

also containing the physical matrix element of operator A4. The two formulae ea. 1451 
and ea. 1501 have a lot of similarities. They are both an expansion in terms of spherical 
functions of different type (Y^ for the canonical and D^, m for the helicity formalism) 
with partial amplitudes which cover the decay matrix element. 

To obtain an intensity it is necessary to include the population of spin states of the 
initial state. For a particle with spin J there ae 2 J + 1 projections to the quantization 
axis, thus leading to a spin-density matrix p of rank 2J + 1 which is usually diagonal and 
has the form 

/ 1 

(52) PMM' = 

V 1 
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With this definition the observed number of events is given by 

(53) /(*)aa'= E A{^ t { V ^)p MMI A{%l{^) 

M,M' ,X C X S > ,A t A t / 

where the summation over A s , X s >, At, Af is implicitly constrained to A = A s — At and 
A' = A s ' — Af 

2'3.3. Relations between Canonical and Helicity Base. Since both approaches deliver 
a complete description of the state, it is fairly easy to go from one representation to 
another. This is important if L and S are good quantum numbers for the system. Then 
the recoupling coefficients lead to symmetry relations among the helicity amplitudes. 

To move from one representation to the other we need the recoupling coefficients. 
They are derived using the properties of the D-functions: 



(UlJMLS\JMX a X t ) = ^|l±i(L0 S(X S - A t )| J(X a - X t ))(sX s t(-X t )\S(X a - A*)) 

With ea . 1541 we obtain a relation between f\ s \ t and a J LS . The recoupling from the canon- 
ical to the helicity representation is 



N jfLx t = ^2(JMX,X t \JMLS)J—(2J + l)(JMLS\M\JM) 



l,s ' P* 



(55) = £ V2L + T{L0 S(X S - X t )\J(X s - X t ))(sX s t{-X t )\S(X s - X t ))a J LS 

L,S 

while from the helicity to the canonical representation is 



a ts 



E (JMLS\JMX s X t )J— (JMX s X t \M\JM) 

A s , A t V Ps 

= (JMLS\JMX s Xt) -j^= fi jXt 



(56) =NjJ2 ^VTTT i 10 5 (A S - A *M A s - A t ))(sA s t(-A t )|S(A s - A t ))// sAt . 

A s , At + 

2'3.4. Symmetry Relations. Reactions involving strong or electroweak interactions 
conserve parity. Parity reverses the direction of r and p while the angular momentum 
remains unchanged. Applying the parity transformation on a single particle state one 
obtains 



(57) 



P\p, jm) — r]\ip + 7r, 7r — jm) 
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in canonical basis and 

(58) P\p, jX) = r,e- m j\ip + tt,it - 0, \p\,j - A) 

in helicity basis with rj being the intrinsic parity (±1) of the particle. For the two-particle 
states the parity transformation leads to the following relations 

(59) P\ JMls) = rnr]2{-l) l \JMls) 
for the canonical and 

(60) | JMX 1 X 2 ) = ]T J^jjji 10 sA| JM)( Sl Ai s 2 (-A 2 )|sA)| JMls) 

1,8 

(61) P|JMAiA 2 ) = r ]1 r l2 (-l) J+s ^\JMls) 



for the helicity basis with r\\ , 772 and 77 being the intrinsic parities of the two daughters 
and the decaying particle respectively. This can be used to derive relations for the helicity 
amplitudes. If parity conservation holds then 

(62) F( iX2 = mm {-l) J+s ^F(_ Xi){ _ X2) 
In the special case, that both particles are identical one obtains 

(63) F( iX2 =%(-!) J F X2 Al . 

2'3.5. Examples. In this section several examples are given to illustrate the actual 
handling of the helicity formalism in the day-to-day business. 

Example: / 2 (1270) — > tttt. The 72(1270) is an isoscalar tensor particle which means that 
the initial state has J G (J re ') = 0+(2++). The two final state pions have 7 G (J R7 ) = 
1~(0 '"). The intrinsic parity of the /2(1270) is even, since 77/ = rj^(—l) 1 and I = 2. The 
total spin s = 2s w is zero. Starting with the definition of the formalism eg. 1501 

KL^NjF^ iX2 Di^,d) 

and using A = Ai — A2 = and J = 2 we get 

(64) A^{^d)=N 2 F^D^{^) 

N 2 F$ = Vb(20 00|00) (00 00|00)a 20 = V5a 20 
1 

(65) A 2 M(<p, t?) = V5a 20 D% (<p, 7?). 
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Remember that this is a short-hand writing for a whole matrix of amplitudes 



(66) 



i4.Q (<p,i?) = V5a 2 o 



The intensity is then derived by the bilinear sum of amplitude and conjugated amplitude, 
weighted by the spin density matrix 



1(0) = ]T A™{^<&) PMM ,Alf*{ V ,$) 



M.M' 



(67) 



P= 5 



Due to cancellation of cost? and sin$ terms, the final result is very simple, it is just a 
constant: 



(68) 



J(0) = \a 2Q \ 2 0^sin 4 i?+ 15sin 2 ??cos 2 t9 + 5 Q cos 2 - ^ 
= | Q.20 1 2 = const 



Example: lo — > jtt. The w is an isoscalar vector particle which means that the initial 
state has 7 G (J K; ') = 0"(1~ ). The final state pion has J G (J re ) = l-(0 _+ ) and the 
photon has I G (J FC ') = 0(l - ~). One should keep in mind that the real photon has no 
longitudinal component. The intrinsic parity of the u> is r/ w = ^^(—l) and the total 
spin s = s w + s-y is 1. Starting with the definition of the formalism eg . 151)1 

= NjFt lX2 D J M \(<p,#) 

and using A = Ai — A 2 = A 7 = Ai and J = 1 we get 

(69) A^(<p,4)=NAD&x(<P,#) 

NxF{ = V3(10 1A|JA) (1A 00|lA)ou = -XJ^au 

-^L. 1 

V2 

(70) Alf&,#) = -xJlanDfafafi). 
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or in matrix representation 

/r/i-D(-i)(^ i-u^ 

(71) Ai ftM?) = -J- -dj ( _ 1} (i?) -<(,?) 

leading to the intensity 

J(0) = Yl A^(^)p MM ,A\^*{^d)5 xxl 
M,M\\.y 

(72) p-i^ 1 i 

which finally collapses again to a constant: 

r . n . 1, l2 / /1-costA 2 /l + cos??\ 2 „sin 2 i?\ 

(73) = i|an| 2 (l + cos 2 z9 + sin 2 tf) = |an| 2 = const 

Example: fo2( an y) ~ * 77- The /o,2 is an isoscalar scalar or tensor particle which means 
that the initial state has I G (J' FC ') = + (0 ++ or 2 ++ ). The two final state gammas have 
I (J 1 * ) — 0(1~ — ). The intrinsic parity of the /2(1270) is even, since i)f = ry 2 (— 1)' and 
I = 2. The total spin s = 2s 7 is two and thus I = 0, 2 are allowed for the /o and / = 0, 2, 4 
are allowed for the f 2 - In the first case we investigate J = 0. Following the definition of 
the formalism ea. l5UI 

a j x % 2 =NjfZ iX2 d j m * x (<p,#) 

and using A = Ai — A2 we get with J = 

(74) 4° A >,tf) = N Fl X2 D° l(v,V) 

NoFg = J2(l0sX\ J\)( Sl \ lS2 (-\ 2 )\ s \)a ls 

l,s 

= vT(0000|00)(lAil(-A 2 )|OA)o o 
+A/5(2020|00)(lA 1 l(-A 2 )|2A)a22 



3«00 + y g«22 

(75) K X2 {^d) = (\^a 00 + Jloaa J D&fol?) = const 
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In the more complicated case for J = we get 



(76) A™[ 2 0) = AT 2 F A ° lA3J D|j A (^ 0) 

(77) N 2 Fl iX2 =^(/0 S A|JA)( Sl A lS2 (-A 2 )| S A)a ;s 

= V5(2000|20)(1A 1 1(-A 2 )|00)a 20 
+V5(202A|2A)(lA 1 l(-A 2 )|2A)a 22 
+\/9(402A|2A)(lA 1 l(-A 2 )|2A)a 42 

(78) = J^a 00 + \f^ a 22 

(79) Af lX2 (<p, 0) = ( J|aoo + J\a 22 ) D° *(<p, 0) const 



Because of the complexity of this result we limit the discussion now to the case where 
J = A = 2. Then a lot of terms disappear and the amplitude contracts to 

A^-i) = ^0s2|22)( Sl ls 2 l|s2K 

(80) = V5(20 22|22) (11 ll|22)a 22 + V9(40 22|22) (11 ll|22)a 42 

/T 1 i 1 

Using the symmetrization because of the two identical particles in the final state one gets 
a fairly simple angular distribution 

(81) A 2M = N 2 (*>!(_!) + F ( _ 1}1 ) DSa^i?) <x £$ 2 (^,0). 

Example: pp -> tttt. The pp(0 _+ ) has I G (J W ) = l+(J rc ) with M = 0,±1. The two 
final state pions have I (J 1 ^ ) = 1~(0 h ). Following the definition of the formalism 
cq.EOl 

and using A = Ai — A 2 = and J = I we get 
o o 

(82) A J 00 M fa 0) = NjF J D{Z fa §) 

NjF °° = V2/TT(/0 00 1 JO) (00 00|00)a /0 = \/2J + la J0 

(83) A 7 M (^, 0) = \727TTaj OJ Dj{;o(<^, - V / 2TTTa J0 < We~ iM ^ 
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The d-functions are not orthogonal, if ip is not observed ambiguities remain in the am- 
plitude and polarization is needed, to disentangle the waves. 

Example: pp — > ojit — > jtt°tt . This is a two step process pp — > ojit and ui — > jn° 
and illustrates how subsequent decays are handled in the helicity formalism. The pp has 
I G (J PC ) = l+iJ^). The intermediate w has J G (J K7 ) = 0~(1 — ), the two pions have 
7 G (J K7 ) = 1~(0~ + ) and the photon has J G (J rc ) = 0(1 ). The angular momentum 
between the recoil pion and the /2(1270) is L = 2. Using the result from eq.[7U|we obtain 

(84) A(f{ 2 (n x , (h) = a(X (o W7r c , n^) = a]^ (n ff o 7 )A{f (o W7r c ) 

= N Ut i F^ D\l x ^ (^o 7 )NppjF Xu0 D^ Xui (f2 W7r o ) 

}1* 



(85) X ^V2m(Z01A u |JA w )a ppiU J 

The interesting feature of this amplitude is, that the helicity constant 0^^0.11 factorizes 
and all helicity amplitudes of the pp system can be modified to 

(86) a p"p,n = a ppji a ^rf),n 

Example: pp — + /2(1270)7r° — * 7r + 7r~7r . This is an two step process pp — * /2(1270)7r° and 
/ 2 (1270) -> vr+vr-. The pp(0"+) has Z G (J- FC ') = 1"(0-+). The intermediate / 2 (1270) 
has ^(J^ 7 ) = 0+(2++) and the two pions have 7 G (J ro ) = 1~((T + ). The angular 
momentum between the recoil pion and the /2(1270) is L — 2. Starting with the definition 
of the formalismea .1501 

K{ 2 ^NjF J XlX2 Di* IX {^) 

and A = Ai — A2 = 0, Jpp = and Jf 27T ° = 2 we get The full amplitude is then given as 
(amplitude tree) 

(87) A^^.na) = A^n^.fW) = ^(^WO^fW) 

= ^pp.O-^OO^OoC^ATrO^^oO-OooC^Tr+TT-P) 

%>,o*oo = ^(20 20|00) (20 00|20)a pPi22 



A Primer on Partial Wave Analysis 



21 




71(20 00|20) (00 00|00)a /2i20 




1 l 



(88) 





with £l xy = (ip, 1?) being the direction of daughter x in the mother system of xy. The 
final intensity turn out to be simple again: 



2'3.6. The pp system. 

Proton antiproton at rest. The proton antiproton system at rest is formed by an incident 
antiproton beam on a hydrogen target. Due to Coulomb scattering the antiproton is 
decelerated according to Bethe-Bloch's formula until it's energy is small enough to be 
caught by a hydrogen atom. The protonium atom (proton antiproton atomic bound 
state) is formed at very high n and I of about 30 where it replaces the electron. The 
system deexcites through slow radiative transitions and through collisions (Auger effect). 
These processes compete with the Stark mixing of the I levels. Due to the fact, that the 
protonium atom is much smaller than a usual hydrogen atom it traverses those and 
feels a very strong electric field which make the levels mix (Day, Snow and Suchcr, 
1960). The average electric field which is seen by the protonium atom varies with the 
density of hydrogen atoms and thus depends on the target density. Since in S-wave 
orbits a annihilation is very probable high density targets (like in liquid hydrogen) have 
a dominant S-wave yield. In low density targets (like in gaseous hydrogen) more P-wave 
is present The advantages of this initial system are manifold 

• varies with the target density, 

• isospin varies with n (deuterium) or p (hydrogen) targets, 

• atomic (incoherent) initial states enable unambiguous partial wave decomposition. 
But there are also some disadvantages like 

• limited phase space and a very small kaon yield. 

The pp system has an invariant mass of mp-p «1877 McV/c 2 . Since a recoil particle 
is needed to ensure energy and momentum balance investigations are only possible up 
to a mass of mueson ~1700 MeV/c 2 . If due to quantum number requirements another 
(heavier) recoil particle is required, the phase space limits arc even more severe. 



(89) 
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In principle the same kind of physics could be investigated using an n beam. This 
has been done by the Obelix collaboration. Since the n as the n is neutral, there is 
now deceleration in material and no neutronium atom. But one can do annihilations in 
flight which is discussed in the next subsection. This is not true for annihilation of an 
antiproton on the neutron of a deuteron, since the deuteron and the antiproton also form 
a bound atomic system. 

The annihilation - being a strong process - conserves all quantum numbers from the 
atomic state. The most important ones are 

• G=(-1) I+L+S , 

• P = 

• C = (-1) L+S and CP = (-l) s+1 . 

The isospin of a pp system can be either or 1. Therefore in the isospin space the 
protonium system mixes with the neutronium and the wave functions for the initial state 
are 

• TfdPP) + l 5n » lms /=0 ' 

• 7f (|PP> - l Hn >) nas I=l -, h=0 and 

• |pn) has 1=1, i3=+l. 

Usually only S- and P-waves are taken into account. In liquid hydrogen the S-wave 
is dominant with a yield of more than 90%. In gaseous hydrogen at NTP the ratio 
between S- and P-wave is around 1. In low pressure targets the P-wave yield may be 
as much as 90% [7]. Since the atomic system lives quite long compared to the strong 
interaction it rotates many times until it finally annihilates. Therefore the helicity of the 
initial antiproton is of no use. The population density is just equally distributed for all 
helicities or in other terms for all values of the magnetic quantum number. 

Proton antiproton in flight. With higher antiproton momenta a higher center of mass 
energy can be reached and thus more massive mesons can be produced. In addition there 
is the possibility to create resonances without recoil particles which makes the analysis 
much easier and less ambiguous, although not all quantum numbers are accessible (only 
fermion anti-fcrmion quantum numbers remain to be producible). The annihilation of 
pp at higher incident momenta is completely different to the situation at rest. There 
are no more atomic levels. As soon as the antiproton is to fast to be captured one deals 
with a scattering process. The number of partial waves rise with increasing momenta 
and detailed investigations show that the maximum needed angular momentum is / ~ 
Pcms 1 200 Me V/c 2 which is almost compatible with statistical models. This leads to a 
huger amount of waves which where most of them are allowed to interfere. But there are 
some constraints which limit the amplitudes to be considered. 
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Since we have a fermion anti-fermion system helicity conservation reduces the initial 
helicities to or 1 with the symmetry relation H„ lV2 = rjj(—l) J H_ U2 , v \ with helicities 
^1,2 for the proton and the antiproton respectively and total angular momentum J. From 
C and CP invariance on gets in addition H\\ = 0ifi + 5+ Jis odd and H_ 1 ^- j — 
if S = and/or J = 0. Since C and CP are conserved interference is only possible 
in between waves where C and CP are equal. This decouples singlet and triplet states 
as well as it decouples even from odd L. In total four incoherent sets of amplitudes 
remain. The final state poses additional constraints so that many amplitudes vanish due 
to other conservation laws. Nevertheless the number of parameters for a partial wave 
analysis is still very high and can easily exceed 100. Therefore it might be possible that 
in many cases the amplitudes have to be simplified with the cost of giving informations 
and constraints away. 

If there is a formation process with two particles in the final state a full featured 
amplitude is usually applicable and produces unambiguous results jS] . But if there are 
three and more particles involved and the number of resonances is large it is necessary 
to reduce the complexity of the fit by integrating out the production process and con- 
centrating on the final state. Many constraints from the initial state are lost that way 
but the number of parameters can be handled easier. The main caveat is the handling of 
coherence. For two fully coherent amplitudes the intensity looks like I = \A + ei<f>B\ 2 = 
| A | 2 + 1 B | 2 + 2 (3? [ AB * } sin ip + Q [AB * } cos ip) . If there is no coherence it is I = | A | 2 + 1 B \ 2 . 
In the case of a mixture of coherent and incoherent terms an effective coherence is needed 
which can be defined like I = \A\ 2 + \B\ 2 +c(3t [AB*] sin (p+% [AB*] cos ip) with c=-2,. . .,2. 
This technique leaves in the possibility to check whether or not coherence between two 
amplitudes is needed at all |||. 

Other anti-nucleon nucleon systems. There are actually two other techniques which have 
been used apart from pp annihilation. One is the nn annihilation by Obelix. This 
reaction cannot occur from an atomic state, but as an in flight reaction. This is done 
by using a antiproton beam to produce antineutrons in the reaction pp — * nn. This 
produces antineutrons in a broad energy range. They annihilate on a deuterium target 
and the neutral events are selected with zero net charge in the final event. The quantum 
numbers can be applied as in the pp case in flight. The usual energy range does not 
exceed 300 MeV/c so angular momenta up to D-w&ve are sufficient to fit the data. The 
other technique is to use a antiproton beam on a deuterium target and by selecting the net 
charge of the event a selection of pp and pn reactions is possible. In addition information 
comes from the spectating nucleon, either a proton or a neutron. The antiproton is in an 
atomic orbit of the deuteron before annihilation, but the levels are somewhat different 
(due to the different reduces mass of the system) and the base system for the waves is 
the center of mass of the deuteron, and not the individual nucleon which takes part of 
the annihilation. In addition, there is Fermi motion of the nucleons inside the deuteron. 
To select the initial angular momentum of the whole system the breakup momentum 
of the spectator is used. If the momentum is small (p spe c <100MeV/c) the process is 
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dominated by S-wave. Due to the momentum range all phase space boundaries (like in 
Dalitz plots) remain fuzzy. 

2'4. Moments Analysis. - The basic idea of this approach is to do a Fourier decom- 
position of of he final state. This is done for each bin of invariant masses of the exclusive 
final state system. This method is explained best by using an example. For an unpo- 
larized target and without measuring the recoil polarization of the neutron in a charge 
exchange reaction of the type 

(90) n~ + p -> AI° + n and M -> a + b 

the differential cross section can be written ^D] as 

where t is the momentum transfer between the incident ir~ and the M° system, M is the 
A/ invariant mass, A p , A„ are the proton and neutron helicities and (p are the polar 
and azimuthal angles of either particle a, or b in the Gottfried-Jackson frame ( 1 ). The 
full hclicity amplitude H\ p \ a (t, M, ip) which is the sum of the terms corresponding to 
all possible intermediate spin j and helicity m states of M° is given by 

oo j 

(92) V^H XpXn (t,M,#,<p) = Y^ £ V2jTT^ p A„,™^oW^- 

j=0 m=-j 

Inserting eq.[521in eq.|^one gets 

M(^) = ^E £ £ V2h + lV2h+Te< m ^- m ^ 

A p ,A n i2,m 2 ji,mi 

The helicity products of relation ea. 1931 can be expressed in terms of the density matrix 

(OA) Jl02 _ _J_ \ " Tjjl* Tjh 

\ I Fmim 2 Z-~t A p ,A„ : mi A p ,A n ,m 2 

A p ,A n 

d 2 a 

where N = — — — . Therefore eq. ESI can be rewritten as 
otoM 

m{ #^ ) = E_ \/2jiTTy277TTp mim2 JU2C l0 W< We< m i-'2)^ 

32,"2jl,™l 

( 1 ) The Gottfried- Jackson frame is the M° frame rotated so that its z-axis is parallel to the 
incident n~ direction. 
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and its integral gives the trace equality 
(96) 



mm 



j,m 



As the strong interaction conserves parity, for the process 1 + 2^3 + 4, the hclicity 
amplitudes satisfy the relation 

(97) #A 3 ,A 4 ,Ai,A 2 = v(- l ) k H-X 3 -\ i -X 1 -X 2 

whre 77 is the intrinsic parity product of the particles 1, 2, 3 and 4 and k — X^( s i + ^0 
is the sum of their spins and helicities. Combining ea. 1971 and ca.1941 one obtains 



(98) /4f m2 = (-l) r 



(— mi)(— m%)' 



In any reference frame where the z-axis is in the production plane, the most general 
angular distribution of the M° system can be written, according to parity conservation, 
as a sum of the real parts 9? [Y^] of the spherical harmonic moments. For simplicity we 
will write Y{ n instead of 3? [Y l m ] 



OO I 



(99) I(t, M, #,V) = NJ2 E WWW V) 



1=0 m=-l 



where (Y™) are the normalized moments. In particular (Y °) = and Y Q ° = 1. 

Using standard relations amongst rotation matrix products, the (Y™) coefficients can 
be expressed as a sum of the real parts of the elements p 3 m{ 2 m2 



j2,m 2 ,ji,m 1 V 

(100) (Im ji(-mi)|j 2 (-m 2 ))(Z0 jiO|j 2 0)5R [p>^ m2 ] . 

Together with eq.[!2]it is easy to see, that the (Y™) can be expressed as bilinear forms 
of hclicity amplitudes. 

Including the factor N in the coefficients (YJ m ),one obtains a new set of coefficients 
tf 1 which depend on the momentum transfer t and the mass M of the M° system. The 
angular distribution of the produced events can then be written as 

oo / 

(101) i prod (t, m, <p) = e E *r*r o*. <p) 

1=0 m=0 



This set of coefficients is obtained by fitting ea. HUll To obtain the different waves, which 
contribute to a specific moment the relation ea. 11001 has to inverted. Since this is a 
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non-linear system of equation there are several solutions and it has to be ensured in the 
analysis process, that the different solutions show the same physical behavior. Otherwise 
the data will not be conclusive. 

An example with only S- and P-wave may serve to illustrate the details. For I < 1 
there are only two components 

(102) = (5)P (cos^) + 3(P)Pi(costf) = (S) + 3(P) cosi? 



3. Dynamical Amplitudes 




Fig. 6. - s-Channel scattering. 



3*1. S- Matrix. - Consider a two-body scattering of the type ab — > cd (like in fig.|SJ|. 
The differential cross section is given in terms of the invariant amplitude M and the 
'scattering amplitude' / through 

< 103 » f? = <sk (I ) ^ = '« n >i 2 

where 'i' and '/' stand for the initial and final states; SI = (6, <fi) denotes the usual 
spherical coordinate system; and s = m 2 is the square of the center-of-mass (CM) energy. 
The qi(qf) is the breakup momentum in the initial(final) system. [The observed cross 
section is in reality the average of the initial spin states and the sum over all final spin 
states — this is suppressed here for simplicity] The scattering amplitude can be expanded 
in terms of the partial-wave amplitudes 

(104) f fi {Q) = i^J+ljT/^JDfc (<M,0) 

qi j 

where A = A a — and \x = X c — \d in terms of the helicities of the particles involved in the 
scattering ab — > cd. Note that this 'scattering amplitude' is a factor of two bigger than 
that with a more common definition (for example, see Section 5.1, Chung |11|). One may 
in addition note that the argument of the /^-function is frequently given as ((f), 6, —(f)) (see 
Jacob and Wick [2] and Martin and Spearman Q2]). Integrating the differential cross 
section over the angles, one finds, for the cross section in the partial wave J, 

(105) aj i =(^](2J+l)\Tf i (s)\ 2 
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Note that T J has no unit; the unit for the cross section is being carried by qf. It is 
necessary to define more precisely the initial and the final states 

(106) |t) = \ab, JMX a X b ) 

(107) |/) = \cd,JMX c Xd) 

where M is the z-component of total spin J in a coordinate system fixed in the overall CM 
frame and the notations {ab} and {cd} designate additional informations needed to fully 
specify the initial and the final states. Because of conservation of angular momentum, an 
initial state in | JM) remains the same in the scattering process. Note the normalization 
(see Section 4.2, Chung |TT] ) 

(108) (f\i) = 5 l3 

In the remainder of this section and in subsequent sections, it is be understood that the 
ket states mentioned always refer to those of ea. 11071 In particular, explicit references 
to the total angular momentum J will be suppressed. Note that, with this convention, 
one has eliminated the necessity of specifying continuum variables such as angles and 
momenta. 

In general, the amplitude that an initial state |z) will be found in the final state |/) is 

(109) S fi = (f\S\i) 

where S is called the scattering operator. One may remove the probability that the initial 
and final states do not interact at all, by defining the transition operator T through 

(110) S = I + 2iT 

where / is the identity operator. The factors 2 and I have been introduced for convenience. 
From conservation of probability, one deduces that the scattering operator S is unitary, 
i.e. 

(111) S = 5 f 5 = I 
32. T- Matrix. - 

3'2.1. Harmonic Oscillator and the Lorentz-Function. 

Harmonic Oscillator. As an academic example toward the discussion about the properties 
of the T-matrix and how to obtain a reasonable parametrization we go back to the 
classical harmonic oscillator, from the Lagrange function 



(112) 



m . o kx 2 

L = -ttX 2 — 

2 2 
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and the Lagrange equations 

d dL dL 

dt dxi dxi 

one can derive the equation of motion, which is in this case the well known formula 
(114) x + w% = 

If the oscillator is damped and periodic forces are applied to the damped harmonic system 
like 



f f 

(115) x + 2\x + ljIx= — cosu R t = — 3R: 

m m 



leads to the time averaged intensity (over a half period) in the proximity of the resonance 

( 116 ) J M = T-l \ 2 m2 - 

4m (uj r - w ) + y 

This function is called Lorentz function and has the phase 5 with 

(117) tan ,5= V™ R „ . 

and has a maximum for = ujq with a value of I = I(uj ) = . 

4m 

Primer in Scattering Theory. In scattering theory a line shape of a resonance is obtained 
by solving the Schrodingcr equation. For two spinless particles with masses mi and m.2 
the interaction can be described by a spherical potential with V(r) = V(\r\) = V(r). 
The Schrodingcr equation is defined like 

(118) - V 2 * + V(r)* = 

with the reduced mass fi = 17111712 _ -phe incoming particle is assumed to be a free 

mi + m 2 

particle, represented by an incoming wave. If the particle moves in z we get 

(119) = e lkz with k= |. 
The incoming wave can be expanded in partial waves 

oo 

(120) *i = £)tf,(r)i5(cos0) 
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where Ui are the solutions of the free (V = 0) radial equation 

1(1+1) 



(121) 



i 2 

dr 2 



^U l (r)+[k 2 -^-^)U l (r)=0. 



The asymptotic behaviors of the wave for kr — > oo is 



(122) 



A:/'- 



,/,:_^ ! (2i + l) 



1=0 



2ikr 



kr — 



In 



In 



Without interaction we get = ^ . Even with interaction the incoming wave remains 
unchanged, while the phase of the outgoing wave is moved by 26 and the amplitude is 
reduced by a factor r/. The asymptotic behavior is 



(123) 



kr- 



t l (2l + l) 



2tkr 

1=0 



i kr 



In 



i 1 kr + 



In 



The difference between the incoming and outgoing wave defines the scattering wave. 



(124) 
(125) 
(126) 



ikr 



1=0 



= ~jT(2Z + l)T,fi(cos0) 



2% 

ikr 



1=0 



„lkr 



1(0)- 



Eq. 11261 contains two important definitions: 



(127) 
(128) 



m e 2l6 i - 1 
2% 



1 OO 

}{0) = T ^(21 + WPl(0OS t 



1=0 



Another important result, is that the total cross section is proportional to the imaginary 
part of the forward scattering amplitude. This is well known as the optical theorem: 



(129) 



^=j9[/(0)] 
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with the scattering amplitude /(#) defined as 

ikr 



*5 = /(#)- 



(130) 



. oo 

^) = -^(2i + l) 



1=0 

oo 



2/ 



: P/(cOS$) 



1=0 



Using = => cos = 1 => Pi (1) = 1 we get 



(131) 



a T = g e + en = 5^(21 + 1)3 [T,] 



To derive explicit forms for cte and 07 we now investigate the imaginary part of the 
T-matrix: 



(132) 



(133) 



3 [T{\ = 3f 



r?/jj - 1 



3 [rji sin(2<5/) + i(— 1)\ cos 2(5/ + 1)] 
(-77/ cos 2(5/ + 2) 

(-77/ COS2(5/ + l+77/ 2 + l-7// 2 ) 



m***' - 1 2 



+ (!-»??) 



From this we get for the elastic cross section (using ea . 11310 



. oo 



(134) 



fe 
4tt 



i=0 
oo 



2i6, i 2 



27 



5>i+i)|t,i 



and correspondingly for the inelastic cross section 



(135) 



1=0 
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For the fully elastic case the T-matrix reduces to the simple form 

e 2lSl - 1 x 
(136) 7] = = e»l sin(fc) 

2i 

sin(^) sin 5; 



e i5i cos (5; — i sin <5/ 

(137) = — ± . 

cot oi — i 

3 '2. 2. Simple and Relativistic Breit-Wigner Function. 



Decay of Unstable States. We assume a state which can be either a resonance in scattering 
theory or a metastable state in an atom. The wave function of the non-stationary state 
of frequency u)r — and lifetime = & can be written as 



i/>(t) = * e _a;fl< e 2r 

r 



(138) 



* e"^e 2 



f as a function of the frequency is derived from a Fourier transformation 

1 



(139) 



V27T 



i){t)e lujt dt 



l\0J-W R + l^-) t 

e V 2 / dt 



(140) 



a; — ujr — i- 



[E R -E)-%\ 



1 i \ w - w fi + 1— J t 



'2tt 



where the constant k is determined by the normalization of the wave function. To get 
the connection to scattering theory, we determine k with the formulae of the previous 
paragraphs. In the case of elastic scattering (the resonance decays only in one elastic 
channel mj + m 2 —> R — > mi + m 2 ) the maximum cross section for a given I is (see 
eq.EH) 



(141) 



(4) max =^(2Z + l). 
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On the other hand the elastic cross section is proportional to the absolute square of the 

4k 2 

wave function oe oc ^"5*. Since ( , I n I'*) max = we get the Breit-Wigner formula with 
r 



(142) V(E) = 



(E R -E)-i± 
The elastic cross section is then 

r 2 

Att — 
(143) a l E = -t-(21 + 1)- 



Breit-Wigner from Phase Movement. In the case of elastic scattering of spinless particles 
via a resonance of spin J = I we start with T = cot g_ t ■ In the proximity of the resonance 

2 



it is well known that S ~ — => cot S w and Er — E(5 — §)■ We expand cot 5 and get 



± cot 5(E) 



E=E T . 



(144) cot 5(E) = cot 5(E R ) + (E - E R ) 

(145) =0 + ( E -E R )(-l) 

where £ was defined as the first derivative of the cot 5. From this we obtain 



(146) T 



cot5-i (E R -E)-i% 
which is only valid for \E — R R \ kT« E r . 

Breit-Wigner from Field Theory. A different way is to start from field theory and to 
to derive a dynamic function, e.g. a Breit-Wigner, from a propagator approach. Let's 

Fig. 7. - Resonant scattering. A pair of mesons undergoes a resonant scattering and appears 
again in the final state. 

suppose we start with a process like the one in fig.0 The transition amplitude, respective 
the T-matrix can be written like 



(147) 
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M, M, M, M, M, 
m 



\m /- m, 



M 2 M 2 + M 2 M 2 M 2 + 
M t M, M, Mi M, 



M, M, M, M 



2 +...-- 



m 

M, M, 



Fig. 8. - T-Matrix from perturbation. 



where the transition is constructed from a propagator and the two couplings at the 
vertices. But that is not the end of the story since more loops can open up like 

1 1 Vi?bVi9 

(148) T = V 12 ——b——V 12 



2 



' E — E E — E (E — E) 

and finally the transition amplitude is constructed from an infinite series of propagations 
(see fig. HI where the strength gets smaller by the factor b for each additional loop: 

(149) T= Vl >V l2 V " W " ■ 



'E -E [E Q — E) 2 {E -Ef 

This can be written as a geometric series 



E - E V E -E (Eo- Ef 

V12V12 ( 1 ^ 
Eo - E \ 1 



E -E . 



(150) Vl2Vl2 



Eo-E-b 
Since b is complex we get 

^12^12 



T = 
(151) 



E a - [b] - E - z3 [b] 
E R -E-i% [b] ' 



Relativistic Breit- Wigner. All discussions so far were based on non-relativistic forms. 
Since most reactions in meson spectroscopy involve relativistic particles a relativistic 
formulae is mandatory for any serious data analysis. Without prove and omitting some 
details here is sketch toward a covariant description. 



34 



K. Peters 



Starting with the optical theorem in a slightly different form 
(152) 3[r( a ,0)] = f ^-a T 



we derive 



(153) T,(s) 



4tt 2 



2 cot 8 — i 

with the inverse T-matrix (dropping I for ease of writing) 

(154) T-\s) = -?=(cot<y- t) = fe(s) - A 

V s V s 

The matrix k(s) will later be introduced as X-matrix (see sec. l3'3.l|) . k(s) can be any real 

valued function. T(s) has a 'peak', e.g. a resonance, if fc(s) = 0. The simplest covariant 

s — s 

function with this property is k(s) = — from which we obtain the relativistic Breit- 

7 

Wigner function 

(155) T(s) 



_ _ 2qj_ 



3'2.3. Centrifugal Barrier. From the covariant description of the decay amplitudes, 
one can show that the mass dependent width of the resonances are proportional to q 2l+1 . 



This is known as the centrifugal barrier. Taking out the phase space factor p — y /, 2q/m 
the remaining width is still proportional to q 21 . This is only valid for low energies (where 
they are very important) and a more general formalism is needed to accommodate for the 
properties of a scattered wave far away from thresholds where the q 2l+1 rule is invalid. 

For fixed-channel orbital angular momentum L these effects are correlated in a simple 
way with the semi-classical impact parameter 

(156) b= [L(L + l)Y/q 

where q is the channel wave number. When b is large (q small), cross sections and partial 
decay widths are relatively suppressed. In any two-body resonance decay channel n 
with orbital momentum L ^ 0, the form of the radial wave function will be determined 
dominantly by the repulsive centrifugal barrier and the channel wave number q n for a 
particle separation r greater than a certain 'interaction radius' (R). The approximate 
wave equation which holds outside R is 

d 2 (b 2 \ 

(157) c^^-ij^p 
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where p = q n r and b n = [L(L + 1)] 2 /q n \s again the impact parameter. The solution for 
eq. 11571 is proportional to the spherical Hankel function 

(158) Ur P % R iCnph^ip) ~C n e 

where C n is a constant. Consider the situation if R is very small compared to b: for 
r < b, the radial probability density grows rapidly with decreasing r. The reason is that 
the outgoing wave is reflected inward by the centrifugal barrier. At r = R, the radial 
probability density is 

(159) [T^{R/b)]- l ^p 2 \hf\p)\\ 

For a decay, we are interested on the lifetime which is involved with the centrifugal 
barrier. In general we obtain a reduced width 7„. In addition the factor q n jm is taken 
out of the physical partial width r„ so that 




(160) 
(161) 



r n = ln ^T?(R/b n ) 



r„(g„) = r 1 



^TMR/bl) 



where m is the mass of the decaying particle and q n /m accounts for the two-body phase 
space. How do they actually look like. The spherical Hankel functions are defined in 
terms of ordinary Hankel functions of half-odd-integer order: 



ni (x) 



(162) 



,(L 2 ), 



7T 2 

2x 

1 

7T 2 



N lH (x) 



For real values of x, h^l^x) — [hf\x)]* . The explicit forms of the first three spherical 
Hankel functions are 

IX 

(163) 4 13 (*) = — 

IX 
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The asymptotic behavior of the spherical Hankel functions for small arguments x <C I 
may be inferred from the forms of the spherical Bessel functions 



(166) ji(x) 

ni(x) lJrl 



(2/ + 1)!! 
-(21- 1)!! 



where the double factorial is defined as a;!! = x(x — 2)(x — 4) The ea . 1163^1651 then 

translate into 



F,(g) = F l (x = ^—)= — 

Qscale ^ \hY\x = 1)| S 



(167) Fj^g) ^ a,e l 



And finally we obtain the damping functions for I = to 3. 

F (x) = 1 
Fi(x) -- 



x+1 



(168) F 2 (a:) 

F 3 (x) 



13x 2 



(x-3) 2 + 9x 



277x 3 



y x(x - 15) 2 + 9(2x - 5) 2 
They are applied to the natural width with the factor Bf which are defined as 

Fl(q) 



(169) Bi{q,q R ) 



Fi{q R ) 



3'2.4. Required Properties for T. The main requirement is the analyticity of T. An 
empirical fact is the appearance of singularities, e.g. poles, which are accommodated 
by allowing a finite number of poles in T, thus T is meromorph (or has maximum an- 
alyticity). In addition T is unitary and fulfills a dispersion relation which is discussed 
below. 



Unitarity. Starting from the optical theorem ea. 11291 it follows that 
(170) %[T l ]=\T l \ 2 
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%[T jj ]=J2\ T v\ 

3=0 



or for multiple channels and dropping L 
(171) 

(172) 

Collecting terms on the left-hand side, we have 



3=0 



(173) 1 - 2iT£ + 2iT 3] + 4 £ I?.T« = 1 - 2*T* )(% + 2 l T* j ) = 1 

3=0 j=0 



This becomes equal to 
(174) 



J=0 



Since the outgoing channels are different they cannot interfere = 0) when i =/= j. 

Therefore, we can write the unitarity condition on the S'-matrix 



(175) 

and obtain 
(176) 



3=0 



Dispersion Relations. The derivation of the dispersion relation requires a lot of complex 
variable analysis, so we present only a sketch how to obtain the relations. A detailed 
description can be found in |13|. 

From the assumption of analyticity we deduce the integral around a closed contour 
(C) for a point inside the contour. This is done with a Cauchy integral 

which can be expressed in terms of the individual parts of the contour as 
Ti(s) 



(178) 



1 

2m 

1 

2m 



SL + ill )-T l (s' -irj) 
as 

-rv\ S S 



ds 



(mi+ni2) 2 



,r { (*' + «7)-T I (« / -M 7 ) i 
1 

S — S 2l7T 



Co J Ci JCi. 
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where Co, C\, C 2 are the infinite circle and the small circles surrounding sl and (mi + 
m 2 ) 2 sl is where the left-hand cuts start and are process dependent, (mi +m2) 2 defines 
the right-hand cuts (see sec. l3'5.l1 for details) We have assumed r\ is small so that it can 
be neglected in the denominators. One can show that if 

|T,(«)|-> |s| — > 00 

(179) < 0((s - sl)- 1 ) s^s l 

< 0((s — (mi + m 2 ) 2 ) -1 ) s — > (mi + m 2 ) 2 

independently of the direction in the complex plane from which the limits are taken, then 
the integrals over Co, Cl, and C2 make a vanishing contribution in the limit as the circles 
get very large or very small respectively; We shall assume this behavior to hold, as it does 
in analytic potential theory Using the hermitian analyticity (e.g. Ti(s-\-ie) — Ti(s — ie) = 
2i3 [Ti(s + is)]) we have 

(180) T l(S ) = ir»,' + ir 

T J-00 S'-S 7Tj(m 1 +m 2 ) 2 S ~ s 

where the imaginary parts are evaluated at s' + in, rj — > (the physical amplitude is 
Ti(s + ie), e — > 0. s is greater than (mi + m 2 ) 2 ) which may be written 

(181) T,(,) = i/ oo^Mf2L^ 

f J-00 s - s-ie 

since S pz(s)] = between two cuts, sl < s < (mi + m 2 ) 2 . 

3'2.5. Factorization. Eq. 11421 and eg . 11461 show that resonances are related to singu- 
larities of the T- matrix. The task of this section is to show that the T- matrix factorizes in 
the case of many open channels. From eg . 11461 follows that the T-matrix has a singularity 
at 



(182) cot 8 = — ^ => E = E R -i- 



2' 



2 

The determinant of the inverse matrix T~ l has a zero at the singularity of T: 



(183) Det 



T 1 ( E = E R -i^ 







After diagonalization of the inverse T-matrix with a unitary matrix U we get 

Ai 

(184) T D l =\ A 2 J Det [Tq 1 ] = AiA 2 

'■ 



A Primer on Partial Wave Analysis 



39 



For a single singularity there only one Eigen- value vanishes. For a singularity of higher 
order several Eigen-values can be zero. For Ai = C(E — Er) we get 



(185) T L 



C(E - E R ) 
A 2 











/ C(E-Er) 



T, 



D 







For Rjj being a matrix with 



R D = lim (E - E R )T D 

E^Ea 



(186) 



= lim 

E^E B 



c 






E-E R 
A 2 





/ 


1 

c 



















V 











the transformation 

(187) R = u T ^U hence % = Y j U ik ^U kj 

° k ° 

proves that R can be factorized: 



(188) R ij =J2^U ik -Lu, 
Therefore also T can be factorized like 

(189) T tJ = v 1 



3'3. K-Matrix Formalism. 



3'3.1. Definitions and Properties. The X-matrix formalism provides an elegant way 
of expressing the unitarity of the S'-matrix for the processes of the type ab — > cd. It has 
been originally introduced by Wigner ^31 an d Wigner and Eisenbud 15 for study of 
resonances in nuclear reactions. The first use in particle physics goes back to an analysis 
of resonance production in Kp scattering by Dalitz and Tuan |16| . A comprehensive 
review is found in 17 . In this paper we give a concise description of the if-matrix 
formalism for ease of reference. Its generalizations to arbitrary production processes are 
covered in some detail. 

The reader is referred to the text book by Martin and Spearman f° r some of the 
material covered in this note. However, one must note that the definitions given in this 
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Fig. 9. - Scattering process using a A"-matrix propagator. The boxes hide the actual process. 



paper arc different from those used by Martin and Spearman. From the unitarity of the 
S, one gets 

(190) T -T ] = 2zT t T = 2z TT* 
In terms of the inverse operators, ea . 11901 can be rewritten 

(191) (T^y 1 -T- 1 = 2il 
which may further be transformed into 

(192) {T- 1 + irf = T- 1 + il . 
One is now ready to introduce the K operator via 

(193) K' 1 =T~ 1 +%I 

From ea. H92l one finds that the K operator is Hermitian, i.e. 

(194) K = K ] . 

From time reversal invariance of S and T it follows, that the K operator must be sym- 
metric, i.e. the if-matrix may be chosen to be real and symmetric. One can eliminate 
the inverse operators in ea. 11931 bv multiplying by K and T from left and right and vice 
versa, to obtain 

(195) T = K + iTK = K + iKT 
which shows that K and T operators commute, i.e. 



(196) 



[K,T} = 
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and that, solving for T, one gets 

(197) T = K(I -iK)~ l = (I -iK)~ x K 
and 

(198) S= (I + iK){I — iK)~ l = (I - iK)~ l {I + iK). 

Note that the T-matrix is complex only through the i that appears in this formula, i.e. 
T _1 has been explicitly broken up into its real and imaginary parts, see equation ea . 11931 
It is also useful to split T into its real and imaginary part. From eg. 11951 one finds, 
noting that if is a real matrix, 

(199) 3? [T] = {I + K 2 'y 1 K = K(I + K 2 )^ 1 

(200) %[T] = {I + K 2 )- 1 K 2 = K 2 {I + K 2 )- 1 

Combining this result with eg . 11901 one finds that the unitarity takes on the simple form 

(201) 3[T]=T*T = TT* 
or, from eg . 11931 one gets 

(202) 3 [T- 1 ] = -I. 

Consider now an isoscalar tttt scattering in S-wave below y/s = 1 GeV. This is a 
single-channel problem and unitarity is rigorously maintained. From eo . 12071 one may 

set 

(203) S = e 2lS 

where S is the familiar phase shift. The transition amplitude T is given, from eo. 11101 

(204) T = e lS sinS 

Note that the factors 2 and i in eo. lllUl make the T attain the simple, familiar form. This 
formula shows that the trajectory of T in the complex plane (Argand diagram) is a circle 
of a unit diameter with its center at (0, i/2). This is the so-called unitarity circle and the 
physically allowed T should remain at or within this circle. The S"-wave cross section is, 
from eg . 11051 

(205) a = (^j sin 2 S. 
The if-matrix for this case is simply 



(206) 



K = tan 6 



42 



K. Peters 



and a pole in K is therefore associated with S — tt/2. 

Consider next a two-channel problem in which the S'-matrix may be expressed as 2 x 2 
matrices. Let Sij be symmetric; then it would take six parameters to represent three 
complex variables: Sn, S12, and S*22- The unitarity relationship 

(207) SikS* jk = Stj 

imposes three independent equations. This shows that the S'-matrix depends on just 
three parameters. It can be shown readily that the matrix elements are 

Su = 77 e Ml 

(208) S 22 = 77 e 2lfe 

S12 = iy/1 - V 2 e 1 * 12 , $12 = ^1 + ^2 



where Si is the phase shift for channel i and 77 is the inelasticity (rj < 1). Note that there 
exists only one inelasticity in the two-channel case. 

Turning to the i^-matrix representation of the T-matrix, let 

« Z 

where K\i = K21 and Kij = real. Then, from eg. 11931 one finds 

(210) T = I ( K " - W K ^ 

y ' 1-D- i(K n + K 22 ) \ K21 K 22 - %D 



where 

(211) D = KnK 22 - K 2 12 . 

3'3.2. Lorentz-invariant Description. The transition amplitudes T as defined in 
eg. II 101 is not Lorentz invariant. The invariant amplitude is defined through two-body 
wave functions for the initial and the final state, and the process of the derivation involves 
proper normalizations for the two-particle states (see Section 5.1, Chung JI]). The re- 
sulting invariant amplitude contains the inverse square-root of the two-body phase space 
elements in the initial and the final states. The Lorentz-invariant amplitude, denoted T, 
is thus given by 

(212) T ij = {p i }if ij {p j }i. 
In matrix notation, one may write 



(213) 



T={pY- T{p}* 
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and 



(214) S = I + 2i {p}i T {pp 

where the phase-space 'matrix' is diagonal by definition, i.e. 



(215) 



and 



(216) 



Pi 



P = 



2qi 

m 



Pi 
92 



and 



P2 



2^2 

m 



The qi is the breakup momentum in channel i. (Here one considers a two-channel problem 
for simplicity without loss of generality.) The unitarity demands that, from ea. 120 II and 
en. 127721 



(217) 
and 

(218) 



Q T = T*pT = TpT* 



It is in this form one encounters most frequently the unitary conditions of the transition 
matrix in the literature. 

The cross section in the Jth partial wave is given by, from ea . 11051 



(219) 



Note that this formula embodies the familiar presence of the flux factor of the initial 
system and the phase-space factor of the final system in the process ab — ► cd. In the 
if-matrix formalism, one allows for p to become imaginary below a given threshold; 
however, the cross section above has no meaning below a threshold, and one could then 
modify the expression above by multiplying it with two step functions: 9{pf ) and 9(pj). 

One may recapitulate the expressions for the differential cross section and its partial- 
wave expansion in terms of the invariant amplitudes T£(s). For the purpose, one defines 
the 'invariant scattering amplitude' 



(220) 



T fi (Q) = ^(2 J + l)T J fi {s)Di;{<t>, 0, 0) 



and the differential cross section is given by 

da fl (A\ ( p f 



(221) 



pi 



l^(0)| 2 
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The initial and final density of states are, with s = to 2 , 



(222) 
(223) 



\ 



m a + m b 



1 - 



m a - m b 



Pf 



m, 



m d 



m 



m c - m d 



m 



in terms of the particle masses involved in the scattering ab 
phase-space factors are normalized such that 



cd. Note that these 



(224) 



1 as to 



The invariant amplitude Tfi(Q) is dimensionless, and has a partial- wave expansion 
ea. 12201 The partial- wave amplitude T^(s) is related to the if -matrix via eg. 12281 and 
unitarity is preserved if the if -matrix is taken to be real and symmetric. It should be 
noted that the formula for the differential cross section eg . 12211 has no 'arbitrary' nu- 
merical factors. The 'conventional' invariant amplitude, introduced in eg . 11031 is given 
by 



(225) 



M fi = l6TrT fi (Q) 



One may consider again the isoscalar ttt: scattering in S-wave below 1.0 GeV. In terms 
of the phase-shift 5, the invariant amplitude is given by, from eq.| 



(226) 



f=-e l5 sm6 
P 



and when substituted into eo. !219l the cross section ea . 12051 results. These expressions are 
very familiar, and they demonstrate clearly the interplay between the phase shifts, the 
invariant amplitudes and the cross sections. 

One can similarly define the invariant analogue of the if -matrix through 

(227) K = {p}i K {p}^ 
translating the inverse if -matrix to (ea . 1193( 1 

(228) R- 1 =f- 1 +ip 
which leads to 



(229) 



T = K + iKpT = K + iTpK 
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and 

(230) TpK = Kpf 
Solving for T, one obtains 

(231) T = K{I-ipK)- 1 = (I-iKp)- 1 K 
and, from ea . 11981 

S=(I + i{p}i K{p}?)(I-i{p}i Kip}?)- 1 

(232) = (I- i{p}i K{p}?)-\l + i{p}^ K{p}i). 

Note that K and p do not commute. The Lorentz-invariant T-matrix is then given by 

1 



(233) T = 




1 - pip 2 D - i(piKn + P2K22) 
where 

(234) D = K n K 22 - Kf 2 . 

3'3.3. Implementation of Dynamics (and/or Resonances). There are two possibil- 
ities for parameterizing resonances in JT-matrix formalism: Resonances can arise from 
constant X-matrix elements with the energy variation supplied by phase space or from 
strongly varying K- matrix elements (pole terms) supplying a phase motion |24|. They 
differ in their dynamical character. The former are assumed to arise from exchange forces 
in the corresponding hadronic channels (molecular resonances), so that great effects are 
expected near corresponding thresholds, whereas the latter (normal resonances) corre- 
spond to dynamical sources at the constituent level, coupling to the observed hadrons 
through decay |24| . The dynamical origin of resonances has to be determined by ex- 
periment. In the approximation of resonance domination for the amplitudes, one has 
therefore 

(235) gtf = ^ £gHggM 

mi, — m l 

R H 

and 

\ - gm{m)g Rj (m) n 



(236) Ka = 



« ~ / , 



where the sum on R goes over the number of poles with masses ran, and the residue 
functions (expressed in units of energy) are given by 



(237) 



9iti( m ) = m R T Ri (m) 
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where gm{m) is real (but it could be negative) above the threshold for channel i. The 
constant A'-matrix elements have to be dimensionless and real to preserve unitarity. The 
corresponding width T R (m) is 

(238) r R (m) = £ 



for each pole R. 

Consider now a normal resonance R coupling to n open two-body channels, i.e. the 
mass m R is above the threshold of all the two-body channels. The partial widths may 
be given an expression 

(239) T Ri (m) = ^= 7 M [B l m (q, qR )] 2 Pi 

and the residue function by 



(240) gm(m) = "f Ri \Jm R T R B l Ri (q,q R )^/p~. 

The B(m)'s are ratios of usual centrifugal barrier factors in terms of the breakup mo- 
mentum in channel i and the resonance breakup momentum q R for the orbital angular 
momentum I: 

(241) B l Ri (q,q R )- h( ' 1 '' 



Fi{q R ) 



Widely used are Blatt-Weisskopf barrier factors [22] ('sec l3'2.3"j) for more details The 7's 
are real constants (but they can be negative) and may be given the normalization 

(242) £7« = L 



In practice, it is probably better to avoid this normalization condition by using the 
parameters 



(243) 9m=lm^rn R T R 

as variables in the fit. The residue function is then given by 

(244) g Ri {m) = g Ri B l Ri (q, q R )yfpi. 

We define a i^-matrix total width T R and the i^T-matrix partial widths T R i by 

(245) f R = J2 r Ri = T R (m R ) = r° ^ 7 % Pi (m R ) 
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From these one finds 
(246) r R = 



(247) t", - r '" 



(248) g° m 



^° R Pi(m R ) ' 

/ m R f Rl 
V Pi{ m R) ' 



Note that the if -matrix width do not need to be identical with the width which is 
observed in an experimental mass distribution nor with the width of the T-matrix pole 
in the complex energy plane. We will refer the former as r o b s , the latter as T po i e . In the 
limit in which the masses of the decay particles can be neglected compared to m R , one 
has T R {m R ) ~ T° R . 

In terms of the 7's and g°'s, the invariant X-matrix now has a simpler form 



(249) fry = 



-y R g R jm R T R B Rt (q, q R )B l Rj (q, q R ) 



r " R 



9m9RjB Ri (g, q R )B Rj (q, q R ) 

< 1 — rr»2 



allowing for the possibility that 7's and g 0, s can be negative. 

Consider now an isovector P-wave 7T7T scattering at or near the g mass. Then the tttt 
elastic scattering amplitude is given by 

(250) K=^%=t^8 

uIq — m 

where mo is the mass of the g and 5 is the usual phase shift. The mass-dependent width 
is given by 

(251) r(m)=f Q^ [S 1 2 

where r is the if -matrix width and q (qg) is the tttt break-up momentum for the irir 
mass m (mo). Neglecting the angular dependence of the amplitude, one obtains 



(252) T = e zS sm6 



m r 



m,Q — m 2 — z?7ior(m) 



[BHq,Qo) 



The first bracket in ea . 12521 contains the usual Breit-Wigner form and the last bracket 
expresses the two-body phase-space factor. In this simple case, the if -matrix width 
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and the observed width are identical. Note that the phase-space factor is absent in the 
Lorentz-invariant amplitude T given by ea. 12121 The q 2 dependence of the amplitude (in 
\B 1 {q^ go)] 2 ) reflects the fact that both the initial and the final hit systems are in P-wave. 
The normalization for the transition amplitude has been chosen such that 

^ +?' 

(253) T = +i and T = — at m = mo. 

P 

It is seen that the invariant amplitude T is not normalized to 1 but to p^ 1 . It is for this 
reason that the Argand diagram is usually plotted with T and not T. 

3'3.4. Some Examples. 



Coupled Channel Case. In this example we investigate the influence of coupling of a reso- 
nance to a second open channel at the T-matrix parameters (phase-shift 5 and inelasticity 
rj). Consider an S'-wave resonance at 1500McV/c 2 decaying into ttit and K K with a total 
width of 100MeV/c 2 . In the first example we haven chosen the i^T-matrix widths to be 
f™ = 20MeV/c 2 and f KR = 80MeV/c 2 ,in the second one we choose T^tt — 80 MeV/c 2 
and T K x = 20MeV/c 2 . The change of the coupling to the channel KK has no influence 
on the inelasticity rj and on the line shape of |T| 2 . The only visible difference is the 
behavior of the phase motion 6 (see fjg. llOf b)). In the case of a strong coupling to the 
KK channel (T K R = 80 MeV/c 2 ) one cannot decide whether there is a resonance or not 
by measuring the first channel only (here 7T7t), if the errors of the phase-shift are large. 
Note the different scales for the phase-shift 6 for the two examples. 



Nearby Resonances. Consider again a tttt scattering at mass to. But suppose there exist 
two resonances with masses m a and nib coupling to the isoscalar D-wave channel. The 
prescription for the if-matrix in this case is that 

/ ,.s T, m a T a {m) m b r b (m) 

(254) K = — -\ j 2 

— m z m£ — m z 

i.e. the poles are summed in the iif-matrix. The mass-dependent widths are given by 

(255) T R (m)=T J^)(±)[B\ q , qR )r 

V to J \qRj 

where r — a or r = b and and T® are the two observed widths in the problem. q r 
is the 7T7T breakup momentum at to = m r . If m a and nib are far apart relative to the 
widths, then K is dominated either by the first or the second resonance depending on 
whether to is near m a or mb- The transition amplitude is then given merely by the sum 

lB 2 (q,q a )} 2 



(256) 



T 



to 2 - im a T a (m) 



m a 

TO 
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-0.6 -0.4 -0.2 0.2 0.4 0.6 



Fig. 10. - 7T7T dominated decays of a resonance at 1500 MeV/c 2 with a total with of 100 MeV/c 2 . 
(light curve KK dominates, black curve nn dominates). In the case of a dominant nn width, 
the phase (a) looks normal in tt-k, while it changes dramatically and gets tiny if the KK channel 
dominates, (b) shows the line shape and (c) the Argand plot for both channels 



[B 2 (q,q b )} 2 . 

In the limit in which the two states have the same mass, i.e. tuq = m a — m b , then the 
transition amplitude becomes 

, 258 n rp = m [T a (m)+r b (m)] 

m\-m 2 - im [T a (m) + T b (m)] ' 

This shows that the result is a single Breit-Wigner form but its total width is now the 
sum of the two individual widths. In case of two nearby resonances eg. 12571 is not strictly 
valid. For a specific example fig. 1111 shows the transition amplitude T from the correct 
equation ea. 12541 and from the approximate ea. 12571 Note that ea. 12571 exceeds the unitary 
circle. 



(257) 



+ 



m b fl 



— m 2 — im b Y b {m) 



(m±) ( jl 

V m J \q b 



Flatte Formula. As a next example we take the isovector S'-wave scattering with the 
ao(980) coupling to the nrj (channel 1) and KK (channel 2) final states. Then the 
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(a) 




1.2 1.4 1.6 1.6 

m [GeV/c 2 ] 



(b) 




(c) 




-0.2 




0.8-0.6-0.4-0.2 -0 0.2 0.4 0.6 0.8 
Re(T) 



Fig. 11. - Nearby resonances in /^-matrix formalism compared to sum of Breit-Wigner. The 
light curves show the sum of Breit-Wigner functions and the black curves demonstrate the 
effect of a unitarized K-matrix approach. The masses are m a =1.275 Me V/c 2 , r°=185 MeV/c 2 
and m a =1.565 MeV/c 2 , r° = 150 MeV/c 2 . Both resonances decay only to -kit. It is clearly visible 
how the sum of Breit-Wigner functions violate unitarity (c), while the K- matrix formalism leads 
to a nice circle. The effect is also visible in the phases (a) which moves backward for while and 
in the intensity (b) where it exceeds unity. 



elements of the invariant fC-matrix ea. 12491 arc 



(259) 
(260) 
(261) 



K u = 



K m = 



7?m r 



7|m r 



7i72TO r 



The normalized couplings are denoted by 7^ and 7I, which are both dimensionless and 
satisfy 



(262) 



7i+72 =1- 
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Then the T-matrix ea .12551 is given as 

m r 



(263) 

If one sets eq. 

(264) 
so that 

(265) 
then 

(266) 



T 



m^-m 2 - im T (puf + p 2 %) \ 7i72 72 



7f 7i72 



9i = 7 



9i + 9i = m or 



T 



9l 
9i92 



9i92 
9l 



Pi9%) 



This is the Platte formula |25|. 

The ao(980) appears as a 'regular' resonance in the irn system (channel 1). The com- 
parable Breit-Wigner denominator, for to near m r , is m 2 . — r 
approximation. One finds, therefore, 



im r T r in the resonance 



(267) 
(268) 



(11 


2 

) 


~\p 2 (m r )\~ 


V7i/ 




pi(m r ) _ 



r 



m pi(m r )"/( 



in terms of the mass m r and width T r . Note that pi's have been evaluated at to = m r 
where T is expected to attain its maximum value. The above formulas give merely a 
good starting point; in practice one must vary Too and To to fit the irq spectrum. The 
ratio (72/71) 2 is an unknown (commonly fixed at the SU(3) value of 1.5), but the shape 
of the square of the amplitudes depends only weakly on this value. Once the ratio is 
fixed, then and 7! are fixed through the normalization condition eg. 12621 Now we 
turn to few examples how if -matrices have actually been used in various analysis. 



Example: irir S-Wave a la Au, Morgan and Pennington. |24j The work from Au, Morgan 
and Pennington combines most of the data on tttt and KK scattering at that time. For 
the I — tttt S'-wave the following parametrization was chosen 

(269) Ka = *-f* V -, ^ r + V cl [-V - l" 

4m|- ^ (s r - s)(s r - s ) ^ 13 [4:m 2 K 

(270) = (s - s a )K t3 
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Fig. 12. - Difference between Flatte and Breit-Wigner. (a) The Breit-Wigner function (grey 
curve) shows a unit circle, while in the Flatte function (black curves) one sees in channel 1 
a drop in elasticity at threshold of channel 2 and a small circle for channel 2 starting at the 
respective threshold, (b) Interestingly the line shape for channel 1 does not change very much, 
while for channel 2 there is a dramatic effect at threshold. 



Alternatively they parametrize also an M-matrix. The M-matrix approach is defined by 
T = (M-ip)- 1 . They use 



(271) 



s - s 



Ei>r,ii>r.j , , 



Am 2 K 



As we will see in most parameterizations of the irir 5-wave a parameter like sq appears 
to satisfy the condition of the Adler zero. 

Example: tttt S- Wave d la Crystal Barrel. This parametrization was used to analyze 
the Crystal Barrel data of the channels pp — > 37r0, 2tt°t] and 2?77r0 with an open KK 




0.7 0.8 0.9 1 1.1 1.2 1.£ 
m [GeV/c 2 ] 



Fig. 13. - Effect of coupling variations in the Flattee formula. The stronger the coupling the 
more dramatic are the effects on the line shape. The solid, dashed and dotted curves represent 
coupling ratios from 71/72=10 to 1 and 0.1. respectively for the two channels. 
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channel being constrained by scattering data [261 12"?] . A if -Matrix formalism was 
applied using 



(272) K tj = E 



9r,igr,jB\{qi)B\{qj) 



with K being a 3x3 matrix with the channels I—tttt, 2—KK and 3=7777. The centrifugal 
barrier was chosen according to the work of Hippel and Quigg [22] ■ For scalars this factor 
is equal to one. 

Example: tttt S- Wave a la Anisovich and Sarantsev. The work of Anisovich and Sarantsev 
is the most sophisticated so far and combines all Crystal Barrel data as well as scattering 
data and other information. The parametrization is very complex and involves a lot of 
parameters to accomodate the different reactions |28| . 

(273) Kij (s) = f £ + U S -±±£o) i^L 

\ ^ mf-s s + s J s + sao 

where is a 5x5 matrix (i, ^'=1,2,3,4,5), with the following notations for meson states: 
l—irir, 2=KK, 3=7777, 4 =7777' and 5= multi-meson states (like dominant four-pion state 
for s <1.6 GeV/c 2 ). The g r ^ is the the coupling constant of the bare state r to a particular 
channel i. The parameters and sq describe a smooth part of the if -matrix elements. 
The factor (s — sa)/{s — sao) is used to suppress the false kinematic singularity at s = 
in the physical region near the tttt threshold. The parameters sa and Sao are kept to be 
small. 




Fig. 14. - Production process. The boxes hide the actual process. The source of production is 
unknown and therefore replaced by a strength parameter. 

3'4. P-Vector Approach. - So far one has considered s-channel resonances, or 'for- 
mation' of resonances, observed in the two-body scattering of the type ab — > cd. The 
if -matrix formalism can be generalized to cover the case of 'production' of resonances 
in more complex reactions. The key assumption is that the two-body system in the hnal 
state is an isolated one and that the two particles do not simultaneously interact with 
the rest of the final state in the production process. 

According to Aitchison j^H], the production amplitude P should be transformed into 
F in the presence of two-body final state interactions, as follows: 



(274) 



F = (I - %KY X P = TK~ X P 



54 



K. Peters 



Or, taking the invariant form, it may be written 

(275) F = (I — ik P Y x P = TK~ X P 

where P characterizes production of a resonance and F is the resulting invariant ampli- 
tude. Note the following relationships: 

(276) F = {p}^ F and P = {p}^ P 

Consider first a single-channel problem, e.g. the isoscalar tttt system in «S-wave below 
K K threshold. Then, the K is simply given by ea. 12061 and one finds 

(277) F — e lS cos 5 P. 

Unitarity demands that in this case the phase of the total amplitude F and the irn 
scattering phase S should be identical. Thus the final-state interaction brings in a factor 
e lS — this is the familiar Watson's theorem [3U| and the production amplitude P has to 
be a real function. It is emphasized that P must have the same poles as those of the 
isf-matrix; otherwise F would vanish at the pole position (S = 7r/2). 

In general, P and F are both column vectors, n-dimensional for an n-channcl problem. 
If the if-matrix is given as a sum of the poles as in ea . 12351 then the corresponding P- 
vector is 

(278) P t = Y &8n{m) 



m D — m 
r R 



and 

(279) Pi = 



E 



0% 9m(m) 



where f3 R (expressed in units of energy), carries the production information of the reso- 
nance R. The constant (3% is in general complex, but it can be set to be real under certain 
assumptions. One should note, that Longacre pHl E2J used complex (3 throughout 
his analysis. This ansatz was heavily criticize by Morgan and Pennington |2(J| 

The P-vector should contain the same set of poles as those found in the i^T-matrix. 
However, according to Aitchison, constant terms (or a polynomial in energy) can always 
be added to the poles in the P-vector (as in the iiT-matrix) elements, without destroying 
the unitarity of the P-vector (or T-matrix), 

(280) P -> P + di 

where the constant di is in general complex. 

In the production process pp — > ir + ir~ir° was described by a production amplitude 
for the p tttt and /2(1270) rnr intermediate states. The data required introduction of 
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an additional constant amplitude which was interpreted as that for direct three-pion 
production. This direct production amplitude can be described in our formalism with 
an additional constant in the P-vector. 

It is often more convenient to rescale /3°'s 



(281) 



f3 R =f3 R ^Jm R T R 



so that /3's are dimensionless. Then the P-vector reads, from eq. 

0mm m R T R B Rl (q, q R ) 



eq. 



(282) 



where, once again, 7's are real but /3's could be complex. If the production process has 

some known dependence on momentum or angular momentum, the production strength 

[3 should be modified accordingly. 

It is instructive to work out the above formula in the case of a single resonance in a 

/3m r 



single channel. Then one has P = 



TYl 



;B l (q,q ) so that, with Kp given by ea. l25UI 



(283) 



F(m) = [3 



m r 



ttiq — m 2 — imoT(m) 



B l (q,q ). 



This is exactly what one writes down for a Breit-Wigner form, except that one has mul- 
tiplied by an arbitrary constant (3 and the centrifugal damping factor B l (q,qo). This 
provides a X-matrix justification of the traditional 'isobar' model. Note that the numer- 
ator is a constant (independent of m) . 

The scattering process may involve particle exchange in the i-channel (see fig. H5(l 
and the left-hand singularities reach into the physical region. Since the formalism does 
not have any left-hand-cuts, they lead to artificial poles in a P-vcctor approach. Since 
these artificial poles depend on the physical processes involved and differ from reaction 
to reaction. If such poles exist, they can be different in the definition of the if-matrix 
and the P-vector. 



Fig. 15. - t-Channel scattering. 



Cahn and Landshoff |18j state that in some approximations the column vector 
(284) Q = R- 1 P and {p}^Q = Q and 



= K~ X P 
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may be considered a constant in a given limited energy range. Then, one has 

(285) F = TQ and F = TQ 

i.e. the two-body final-state interaction may be expressed as a product of the T-matrix 
and a constant column vector. The Q-vector is devoid of the threshold singularities (i.e., 
no dependence on p) and does not contain pole terms. It therefore depends in general on 
s = m 2 only. In a single-channel problem, e.g. the isoscalar irir system in S'-wave below 
1 GeV, we now obtain, instead of ea. 12771 derived for the P vector approach: 

(286) F = -e lS sm5 Q 

P 

This amplitude contains the familiar scattering amplitude e lS sin 5. 

The P- and Q-vector approaches, even though both being approximations for produc- 
tion of multi particle final states, correspond to different interpretations of the physical 
processes. For clarity we consider a specific reaction, e.g. channel 1: pp — > ttttt] and 
channel 2: pp — > KKr\ from which we want to extract information on (jnt) interactions. 
In the Q- vector approach, the amplitude is given by F± — TnQ\ 4- TyiQ2- the ttit 
system is produced with an amplitude Qi(s). The two-pion interactions are then taken 
into account by the scattering amplitude T\\. Alternatively, a K and a K are produced 
with amplitudes Q2 which scatter via T12 into the outgoing pions. This picture has to 
be contrasted to that which one may have in mind in applying the P-vector approach, 
F = (I — iK)~ 1 P. A resonance is produced with an amplitude P, the term (I — iK)~ x 
may be considered as propagator for this particle which then decays. 

3'5. N/D-Method by Example. - The parametrization in terms of if -matrices with or 
without production vectors is only one approach to describe two-body scattering. Under 
the assumption of unitarity and maximum analyticity it is possible to derive the N/D 
method from the dispersion relations. The amplitude T; for a given angular momentum 
fulfills a dispersion relation derived from the optical theorem. This is an integral equation 
for which has no explicit solution. But T; can be expressed in terms of a nominator Ni 
and a denominator Di 

(287) Tt(a) = 

Di(s) 

where Ni contains only left-side singularities and Di only right-side singularities. The 
functions A; and Di are correlated by a system of integral equations. The advantage 
is, that this system can be reduced to a differential equation of the the Fredholm type, 
which is known to be solvable. The solution is a infinite series, where the actual form 
depend on the integration kernel. In practice this method is realized by a finite sum, 
with parameters being tuned to fit the experimental data. 

The work by Anisovich, Bugg, Sarantsev and Zou |35| may serve as an example for 
such an analysis. In order to fit the Crystal barrel data on pp — > 3ir°, and 27777° they 
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used the following Ansatz below 1.1 GeV/c 2 

Kij = ( i^i) M ™ +a ij+ bij s 

V s J \sa - s s B - s s c - s 

(288) N^(s) = N n (a) = (ci + c 2 s)K n + ip 2 {c 3 + c 4 s)(X 11 X 22 - K 12 K 21 ) 
N vv (s) = N 22 (s) = Cl K 22 + ip 2 c 3 (K n K 22 - K 12 K 21 ) 

where pi are the usual phase space factors and Oy, bij and a are arbitrary complex 
parameters and s» are pole positions, a,, . . ., 7j are coupling constants. 

3"5.1. What is a Resonance. So far we collected several approaches to obtain a 
dynamical function. But all formalisms are only parameterizations which are more or 
less motivated by models. The main physical content is hidden in the T-matrix, however 
it was obtained. 

The instructive Breit-Wigner function is well suited for low statistics, where compli- 
cated methods fail because of the many unconstrained parameters, or in the case of very 
well separated resonances. In this case the fitted mass and width are easily mapped to 
the physical properties of the state. In general this is more complicated and lead s to 
the investigation of the singularities of the meromorph T-matrix. 

The X-matrix poles as well as other parameterization deliver only in rare cases sin- 
gularities which are near to them of the T-matrix. This is for example the case for 

• resonances which are very far apart, 

• minimal coupling to sub-dominant channels and 

• far away from thresholds. 

If at least one of those is violated interference terms appear which move the unphysical 
poles. 



Sheet 3[<2i] Q[q 2 ] 



I + + 

II - + 
III 

IV + 



Table II. - Definition of Riemann sheets (2 channels). 



The energy of a particle is function of p 2 . The inverse function, the root, has no 
unique solution. The sign can be positive or negative. If the momentum is a complex 
number, the two roots lie tt apart on the same circle. This is easily illustrated for w — z 2 . 
lfz = e l V then w = e 2 ^. If ip moves only from to n then w moves around the whole 
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circle. If if moves on w starts a second turn. The difference between the lower and upper 
function value is 



(289) 



= ± 



/ v 


-*-\ 


in 

— cosh — 




he 2 


V 




2 



1(3=0 



^0. 



To resolve this ambiguity one introduces for each threshold involved in the process a 
new pair of Riemann sheets, which represent each one possible root of the breakup 
momentum. The convention is (1) that every crossing of the real energy axis involves a 
sign change for the open channel and (2) that for every passed threshold the sign of the 
threshold channel changes. Please see tbl.lTTl and fig.EDfor illustration. This means for 
a two channel problem 

(290) (r 11 ^ 1 = {r Ii y 1 + i P2 

where (III) and (II) denote Riemann sheets according to the definitions in tbl.[n] The 

Im(E) 



I 



IV 



(a) 



II 



Re(E) 



III 



KK-threshold 
Tut-threshold 
Im(q 2 ) 



II 



(b) 



KK-threshol 




I 



m "*Re(q 2 ) 



Fig. 16. - Numbering of Riemann-sheets (2 channels). 



different sheets in the complex energy and momentum plane are shown in fig. El The 
different sheets are separated by cuts. The convention is that the cuts lie on the real 
axis and start at the threshold. In s-channel it is common to defines the real axis 
above the threshold to be the right-hand cut (RHC). It contains all singularities of the s- 
channel scattering processes. The t- and u-channel can have resonances and bound-states 
too, which depend upon the reaction under investigation and are defined as left-hand 
cuts (LHC), since they appear usually below threshold (see fig.E|f or illustration). One 
should note, that these singularities are not taken into account in i^-matrix formalisms. 
Furthermore, the LHC may reach far out into the physical region, where they might 
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left-hand cuts 



4m -4m; 



4m 2 -nn i 2 



Im(E) 



Re(E) 



right-hand cuts 



Fig. 17. - Illustration of left-hand and right-hand cuts. 



amplitude. To clearly identify resonances it is mandatory to investigate the convergence 
radius of the pole. The convergence radius has to reach out to the real axis. Unfortunately 
the amplitude on the real axis is dominated by the nearest singularity and effects from 
the LHC if they influence the physical region. In summary this means, that the question 
of resonances is translated into the question of how many poles does the T-matrix have 
in the complex energy plane and what are their properties. 

The basic idea is simple. The T-matrix contains a complex denominator D. A 
singularity appears if D vanishes. The mass of a resonance is then the real part of this 
point and the width is twice the imaginary part 



(291) 



which has a similar singularity as a non-relativistic Breit-Wigner. The fundamental 
statement is fact that the relation between the position of the pole and the according 
Riemann sheet characterize the resonance [201 ■ The cross section (real axis) is mainly 
influenced by the nearest pole (see fig. ll8fl . In principle all poles influence the real axis, 
which is easily seen if they are close together. In general the following approximation 
holds near the threshold (at the boundary between sheets II and III) 



(292) 



1 



^iii 



where T r is usually larger than . The empirical principle of meromorph functions 
for T (maximum analyticity or analytic with a finite number of poles) requires that only 
physical singularities are allowed in the analytic continuation into the complex energy 
and momentum plane like 



• kinematic singularities which are defined via unitarity relations. They appear at 
every threshold, where the two-body breakup gets real, 
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• dynamic singularities due to the interaction, like in exchange processes, and 

• poles which correspond to resonances or bound states. 




Re(q) 



Fig. 18. - From a pole to a real Breit-Wigner. The nearest pole dominates the properties of a 
point on the real axis. 



The time dependence of a decay amplitude is the Fourier transform of an exponential 
decay (oc e^ r ). Thus, the imaginary part of q r is negative, which means, that the pole 
lies in the unphysical half of the momentum plane (3 [q] < 0, sec fig.lliJfa'). In contrast to 
that, a bound state lies on the positive imaginary axis. Due to unitarity and hermitian 
analyticity the analytic continuation leads to 



(293) 



T l (q)=T l *(-q*) 



and 



T l (s)=T l *(s*). 



Therefore we have two poles for a resonance symmetric around the imaginary axis. In 
the energy plane they lie symmetric around the real axis. 

I m ( E ) physical plane - Sheet I 



1 1 



Re(E) 



lm(q) 



(a) 



Im(E) unphysical plane - Sheet II 
3* 

— • — ► 
2 #3 



Re(E) 



(b) 



2 

O 

3 3 



*Re(q) 



Fig. 19. - Singularities in the (a) complex energy and (b) momentum plane. (1) bound states, 
(2) anti-bound states and (3) resonances. 



Let's go back to poles at threshold, since they are the most complicated ones and 
the respective resonances remain a mystery in meson spectroscopy. If the poles are far 
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away from thresholds, they are no really affected by them (see fig.l20b). This is true 
even for the p(770) which lies exactly between two threshold. For the Oo/ fo at the KK 
threshold, the situation is more puzzling. Without thresholds, the poles on sheet II and 
III are usually identical and lie on top of each other. At threshold, they do not and it is 
important to verify that the two poles which one finds near the threshold correspond to 
each other. A resonance can only be quoted if the poles belong to each other. Since the 
two sides of the threshold correspond to a different sign for p (sheet II: (iyfp) 2 = —p and 
sheet III: y/p — +p) it is sufficient to show, that while changing p smoothly to — p the 
poles can be connected (see fig.l2Ub). 



p(770) 



IV 



II 



f (980) 



III 

IV 



II 
I 



III 
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IV 



(a) II III 

Im(E) f Tut-threshold 

KK-threshold 




Re(E) 



(b) II III 

Fig. 20. - Normal poles (a,c) and shadow poles at threshold (b). 
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